Spin liquid phase in a S = 1/2 quantum magnet on the kagome lattice 
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We study a model of hard-core bosons with short-range repulsive interactions at half-filling on 
the kagome lattice. Using quantum Monte Carlo (QMC) numerics, we find that this model shows 
a continuous superfluid-insulator quantum phase transition (QPT), with exponents z = 1 and 
v ~ 0.67(5). The insulator, I*, exhibits short-ranged density and bond correlations, topological 
order and exponentially decaying spatial vison correlations, all of which point to a Z2 fractionalized 
phase. We estimate the vison gap in I* from the temperature dependence of the energy. Our results, 
together with the equivalence between hard-core bosons and 5 = 1/2 spins, provide compelling 
evidence for a spin-liquid phase in an easy-axis spin-1/2 model with no special conservation laws. 
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Spin liquid states, or strongly correlated quantum 
paramagnets that preserve all lattice symmetries, were 
proposed long ago as plausible candidates for the ground 
state of frustrated quantum magnets In recent years, 
several frustrated antiferromagnets have been discovered 
which appear to either have a spin-liquid ground state 
(such as the S = 1 kagome lattice system NdaGasSiOu 
[2| or the S = 1/2 triangular lattice Mott insulator k- 
(ET) 2 Cu2(CN) 3 [§]) or are proximate to a spin-liquid 
phase (such as the S = 1/2 stacked triangular magnet 
CS2CUCI4 [H), thus lending support to this paradigm. 
On the theoretical front, there has been considerable 
progress in understanding the effective field theories and 
properties of such spin-liquid phases [B|, @, 0, HI , showing 
that the excitations in this phase carry fractional quan- 
tum numbers and interact with emergent gauge fields. 
However, there is the pressing issue that most micro- 
scopic models which can be shown to have a spin-liquid 
phase have either very unusual interactions [9j or lo- 
cal Hilbert space constraints. An example of the lat- 
ter are particular triangular [13] and kagome lattice [HI] 
quantum dimer models which have dimer-liquid ground 
states but are not directly related to any microscopic spin 
Hamiltonians. 

In a significant development, Balents et al 12] pro- 
posed a model of S — 1/2 spins on the kagome lattice 
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where Sq = J2i<=o * s a sum over ^ e slx s P ms on 
each hexagon of the kagome lattice unit cell, J2q de- 
notes a sum over all hexagons on the lattice. This model 
is easily seen to be a short-range anisotropic XXZ model, 
with only the first, second and third neighbor interac- 
tions being nonzero and equal to each other. Analyz- 
ing this model [3 for Jj_ < and J Z /\J±_\ > 1, they 
showed that it directly maps onto an effective triangular 
lattice quantum dimer model with three dimcrs touching 
each site. In spin language, this effective model takes 
the form of a ring-exchange model, with an exchange 



scale Jring = J\/J z , which describes quantum dynam- 
ics in the Hilbert space with Sq = on each hexagon, 
the local constraint arising from energetic considerations 
at large J z /Jj_. Supplementing this model with an ad- 
ditional four-site (Rokhsar-Kivelson (RK) [13]) potential 
term of strength U4 was shown to lead to a spin-liquid 
for V4 = J r ing, which was argued to be stable for small 
deviations va < J r i ng . Later exact diagonalization (ED) 
numerics [l4j showed that the ring-exchange model ap- 
pears to be in this spin-liquid phase down to V4 = 0, but 
only system sizes upto 20 unit cells could be explored. A 
numerical study of a related rotor model on a decorated 
square lattice [l5( also led to a Z 2 spin liquid, but the 
corresponding S = 1/2 model was not studied. 

In this paper, we show the presence of a spin-liquid 
phase in a S = 1/2 model on the kagome lattice, by 
studying numerically the Hamiltonian {T]) with J± > 
0, J z > 0. Since the ring-exchange model is independent 
of the sign of J± , we expect to recover the physics of the 
ring-exchange model at large J z , but without imposing 
Hilbert space restrictions by hand. On the technical side, 
the choice of J± > eliminates the QMC sign problem. 
This permits us to study very large systems at very low 
temperature using a generalized stochastic series expan- 
sion QMC algorithm 3, 17] and, thus, to go significantly 
beyond earlier work on this class of models. 

By the mapping between 5=1/2 spins and hard- 
core bosons, model ([T]) is equivalent to a hard-core boson 
model with short-range repulsion at half-filling 
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where b\(bj) is the boson creation(annihilation) operator, 
t = J± > is the hopping amplitude, V = J z > is the 
repulsion strength, m = b\bi is the number operator, and 
fi = 12J Z is the chemical potential. The hopping term 
connects only the first, second and third neighbors. We 
begin by showing that model ^ exhibits a superfluid- 
insulator transition at (V/t) c « 19.8. We then turn to 
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FIG. 1: (color online). The superfluid density p s versus V/t, 
at T = t/Q for different linear system sizes L. 



the nature of the insulating phase and show that it is 
a topologically ordered Zi Mott insulator. We thus ar- 
rive at the result that for Jj_ > 0, model {]]) exhibits a 
ferromagnetic phase at small J z and a Z-i fractionalized 
spin- liquid at large J z , separated by a continuous QPT 
at {J z /J±) c « 19.8. 

Superfluid-insulator transition. — For small values of 
V/t, we expect the ground state of @ to be a super- 
fluid {ST). We confirm this by measuring the super- 
fluid density p s through winding number fluctuations 
W ai 2 in each of the two lattice directions, with 
p s = (l/2/3t)((PF Q 2 i ) + (IF Q 2 2 )), where /3 is the inverse 
temperature. As seen from Fig. [1] for small V/t, p s is 
large (its value agrees with mean field estimates [19|). p s 
decreases with increasing V/t, eventually vanishing for 
V/t > 20 suggesting a phase transition to an insulat- 
ing phase (I*)- This behavior is in sharp contrast to a 
nearly identical model where the hopping and repulsive 
interactions are restricted to the nearest neighbor only 

in that case a uniform superfluid persists for arbitrar- 
ily large V/t [2fjj . The absence of a jump in p s in Fig. [T] 
suggests that the ST — I* transition is continuous. 

In the vicinity of a continuous QPT, we expect p s = 
L- z F Ps {L x l v (K c - K),13/L Z ), where F Ps is the scaling 
function, L is the linear system size, z the dynamical 
critical exponent, v the correlation length exponent, and 
(K c — K) oc (V c — V)/t is the distance to the critical 
point. Thus if we plot p s L z as a function of V/t at fixed 
aspect ratio f3/L z , the curves for different system sizes 
should intersect at the critical point. The inset of Fig. [2] 
shows such a plot for z = 1 and (3/L — 3/(4i) with 
a clear crossing point at (V/t) c w 19.8. To obtain the 
correlation length exponent v, we plot p s L as a func- 
tion of [(V/t) c - V/t]lM v for different system sizes. It 
follows from the above scaling relation that the curves 
for different system sizes should collapse onto a universal 
curve F Ps for a properly chosen (V/t) c and v. In Fig. [21 
we show such a data collapse for (V/t) c = 19.80(2) and 
v = 0.67(5). The error bars are estimated from the sta- 
bility of the collapse towards varying the parameters. To 
summarize, we find a continuous ST ~X* transition with 
exponents z — 1 and v ~ 0.67(5). We next examine the 
nature of the insulator X* . 
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FIG. 2: (color online). Data collapse of the superfluid density 
p s for /3/L = 4/(3t), (V/t) c = 19.80(2), and u = 0.67(5). 
Inset: Finite size scaling of p 3 for /3/L = 4/(3t). 



X* does not break lattice symmetries. — For a system of 

l22fl show that an 



bosons, flux-threading arguments [21 



insulating state at half-filling could either be a conven- 
tional state with broken lattice symmetries or must nec- 
essarily have topological order (which means the ground 
state degeneracy depends on the topology of the system). 
We rule out the first possibility here by studying corre- 
lation functions in the insulating state. We look for sig- 
natures of diagonal (density) ordering by studying the 
equal-time density structure factor S(q)/N = (Pq T Pq T ), 
where p qr = {1/N) ^ pi T exp(iqr-j) and pi T is the bo- 
son density at site i and imaginary time r. Fig. [3] com- 
pares a contour plot of the equal time structure factor 
to that of the classical model (with J± = 0) which is 
known to have short-range correlations [23j . The striking 
similarity between the two suggests short-ranged density 
correlations in the full quantum ground state of the in- 
sulator. We confirm this via a careful examination of 
the equal-time structure factor on large system sizes. 
Upto L — 48 we do not find any Bragg peaks which 
rules out the possibility of density ordering even with 
moderately large unit cells. We also examine for possi- 
ble bond ordering in X* by computing the bond-bond 
structure factor Sb(q)/N = (B^B^), where £? qT = 
(1/N) J2 a Bar exp(iqr Q ) is summed over the bond index 
a connecting spins i and j, and B a ^ j^ T = t{b\bj + bib^) 
is the off-diagonal bond operator. Again, we find no sig- 
natures of ordering from the bond-bond structure factor. 

Finally, the ST — T* transition appears to be contin- 
uous, rather than first-order which would be expected 
if we had a conventional transition between a superfluid 
and a broken symmetry insulator. This argument, to- 
gether with the correlation function results, strongly fa- 
vor a non-broken-symmetry ground state for X* . There- 
fore, we next examine the second possibility, that X* is a 
featureless topologically ordered Mott insulator (a "spin- 
liquid" in magnetic language). 

Topological order in X* . — On a lattice with periodic 
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FIG. 3: Contour plots of the equal-time density structure 
factor (L = 24) for the classical model (with t=0) [23| (left 
panel) at T — and in the quantum insulator (right panel) 
for V/t = 20, T = i/12. Axes range from -2tt to 2tt. 



boundary conditions in both lattice directions, the sub- 
space of configurations with no = 3 on every kagome 
hexagon, i.e. the "dimcr subspace" which dominates the 
boson wavefunction for V — > oo, has topological sec- 
tors defined by having, for each lattice direction 01,2, 
an odd (or even) number of bosons on each row ( "parity 
sectors"). In this torus geometry and in the restricted 
Hilbert space, we thus find four topological sectors la- 
belled by row and column parity, such that any local 
move which obeys the local constraint no = 3 leaves the 
sector unchanged. In order to change from one sector to 
another one needs to make highly nonlocal boson moves 
over loops which wind around the lattice. 

The row/column parities are however not well-defined 
in model @ since at any finite V, no matter how large, 
there will be a small density of hexagons with no =/= 3 
24j . We have checked in our numerics, where we start 
from a configuration in the dimer subspace with a certain 
row/column parity, that the QMC algorithm generates a 
small density of particle-hole pair defects in equilibrium, 
so that the quantum ground state no longer lies in the 
"dimcr subspace" . However, for a large enough system 
size at a given value of V/t (L > 10 at V/t — 26), we 
find that our simulations with the longest accessible MC 
steps do not lead to any non-local boson moves which 
wind around the lattice. Thus the winding number iden- 
tically vanishes, and each configuration in the simulation 
which lies in the dimcr subspace belongs to the same 
parity sector as the initial configuration. The full ground 
state accessed by the QMC is, in this sense, perturba- 
tively connected to the initial parity sector. This means 
that four topological sectors continue to exist, and we 
can label them simply by the row/column parity of that 
component of the ground state wavefunction which lies 
in the dimer subspace. 

For a topologically ordered insulator, the ground state 
energy should be identical in each of the four topological 
sectors (on a torus) leading to a ground state degener- 
acy of four. We have computed the energy of the four 
ground states by starting our simulations from configu- 
rations in the dimer subspace lying in four different par- 




FIG. 4: (color online). Equal-time vison-vison correlation 
function for L=24 and T=t/18 in the insulator at V/t=20.5, 
showing exponential decay with a length scale £=1.43(5). 



ity sectors. We find that they are equal within statistical 
errors, which is strong evidence for topological order [25j |. 

Spatial vison correlations in I* . — A second signature 
of Z-2 fractionalization is that visons, which are gapped 
Z<2 vortices in the effective field theory description, should 
have exponentially decaying spatial correlations. The 
spatial vison-vison correlation function is the expecta- 
tion value of a string operator in terms of the spins. For 
the ring-exchange model (valid for V/t — * 00), it takes 
the form C vv (r y ) = |(0| Uk=i e"" fc |0)|, where |0) 

denotes the ground state and the product is along some 



path on the kagome lattice that contains an even num- 
ber of sites, starts at site i, and ends at site j, making 
only "±60°" turns to the left or right, rik = 0, 1 is the 
number of bosons at site k. The absolute value of the 
product is path-independent in the dimer subspace, and 
it is expected to decay exponentially in the topologically 
ordered phase. In model |T]) at finite V/t, ground state 
no longer lies entirely in the dimer subspace, but will 
mix in configurations with particlc-hole defects. Corre- 
spondingly, the vison operator must include correction 
terms to the above definition since it otherwise becomes 
path-dependent in the presence of such defects. Both, the 
corrected ground state wavefunction and the vison opera- 
tor, may be determined order-by-order in a perturbation 
expansion [l9| . undoing the unitary transformation that 
leads from model ([2]) to a ring-exchange Hamiltonian. 

Here evaluate the leading term, correct to 0(t/V), of 
C vv (rij ) by computing the above string expectation value 
in the true ground state projected into the dimer subspace 
[lfH . As seen from Fig. [4j C vv (rij) decays exponentially 
in I* , again signaling topological order in I* . At V/t — 
20.5, we estimate a "decay length", £ = 1.43(5), which 
is comparable to its value at the RK point of the ring- 
exchange model [HI, £ p» 1, and to that found by ED 
[3], £ w 1.7, in the ring-exchange model with V4 = 0. 

Vison gap. — To provide further evidence for gapped 
vison excitations, we display the temperature depen- 
dence of the system energy per site and compressibility 

K = \X£2i n i) 2 ^ m Fig. 03 The energy exhibits a two- 
step decrease upon lowering temperature, with a distinct 
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our model is not of direct relevance to known Kagome 
magnets, our work represents significant progress towards 
understanding realistic Hamiltonians by showing a spin 
liquid phase in a model with only local interactions and 
no special conserved quantities. 
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FIG. 5: (color online). Energy per site E and compressibility 
k versus temperature for L = 24 (V/£=20.5). The energy rises 
exponentially at very low T (see inset) with a wide intermedi- 
ate plateau from T/t ~ 0.5 — 3 (see text). The compressibility 
is zero (within error bar) at low T, rising only at T ~ 3t once 
gapped charge (spinon) excitations become relevant. 

intermediate plateau, before vanishing exponentially at 
a very low temperature. We identify the first drop in 
energy with a freezing out of charge fluctuations below 
a charge gap scale, which we confirm by the sharp de- 
crease in k at this temperature (also shown in Fig. [5]). 
The plateau then corresponds to a regime where the sys- 
tem dominantly explores configurations with no = 3 on 
each hexagon. (The large difference between the total en- 
ergy of the ground state and that in the plateau regime 
can be explained by the macroscopic entropy density of 
multiple vison excitations [l9(.) Finally at the lowest 
temperature, the system begins evolving into the spin- 
liquid ground state. Heating up from T = 0, we therefore 
identify the lowest energy excitations out of the ground 
state as vison-pair excitations (since the charge gap is 
much larger). The temperature dependence of the en- 
ergy then leads to a rough estimate of the single vison 
gap; for V/t = 20.5, we find E v /t ~ 0.35(15). 

Conclusions. — We have studied the T = phase dia- 
gram of a hard-core boson model with short-range repul- 
sion on the kagome lattice using QMC. We find a con- 
tinuous QPT from a superfluid phase to a featureless 
Mott insulator with increasing repulsion. In magnetic 
language, this corresponds to a ferromagnet-paramagnet 
QPT with the paramagnet being a correlated quantum 
spin-liquid. This spin-liquid supports gapped visons and 
exhibits topological order characteristic of a Z2 fraction- 
alized phase. The apparently 3D XY exponents at the 
T* —ST transition, z = 1 and v — 0.67(5), are consistent 
with the QPT arising from condensation of fractionalizcd 
charge- 1/2 bosons [15(. We have estimated the vison gap 
in the spin liquid from the temperature dependence of 
the energy. Finally, the spin liquid phase is fully gapped 
and thus stable to weak perturbations away from the spe- 
cial choice of exchange couplings in our model. Although 
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